Small-scale structure formation properties of chilled sterile neutrinos as dark matter. 
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We calculate the free-streaming length and the phase space density of dark-matter sterile neutrinos 
produced from decays, at the electroweak scale, of a gauge singlet in the Higgs sector. These 
quantities, which depend on the dark-matter production mechanism, are relevant to the study of 
small-scale structure formation and may be used to constrain or rule out dark-matter candidates. 
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I. INTRODUCTION 

The sterile neutrino is an appealing dark-matter candi- 
date. The gauge singlet fcrmions are universally used in 
models of neutrino masscs^l]; the seesaw Lagrangian [2] 
built out of gauge singlet and non-singlet neutrinos can 
explain the observed neutrino masses and mixings for a 
wide range of Majorana masses. If some of these masses 
are small, the corresponding new degrees of freedom ap- 
pear in the low-energy effective theory as sterile neutri- 
nos. Sterile neutrinos with masses in the keV range can 
account for cosmological dark matter [1, 0, S S SJll 
and can explain the observed velocities of pulsars [9|. 
They can also play a key role in baryogenesis and in 
the formation of the first stars [ll| . Unlike many other 
candidates for dark matter, the sterile neutrinos have a 
nonzero free-streaming length that depends on their mass 
and the production history 0, H, HI] . The non-negligible 
free streaming has an observable effect on the structure 
formation. In this paper we will quantify the related 
properties of dark matter in the form of sterile neutrinos 
produced at the electroweak scale. 

Dark matter with a negligibly small free-streaming 
length, called cold dark matter (CDM), is consistent 
with the observations of large-scale structure and mi- 
crowave anisotropy. Warm dark matter (WDM), with a 
free-streaming length smaller than 1 Mpc, fits the large- 
scale structure equally well. The difference arises on the 
smaller scales. CDM predictions for the structure on 
small scales have been studied numerically, and a num- 
ber of inconsistencies have been reported between the 
predictions of CDM and the observations [H, [13, 
In contrast, WDM solves these problems by suppressing 
the structure on the small scales. It is possible that these 
discrepancies could go away as the simulations and the 
observations further improve. However, it is also possible 
that we are seeing the hints of dark matter in the form 
of sterile neutrinos. 

In this paper, we focus on the small-scale structure for- 
mation properties of dark-matter sterile neutrinos, pro- 
duced, at the electroweak scale, from decays of a gauge 
singlet in the Higgs sector. The Higgs singlet vacuum ex- 
pectation value (VEV) gives rise to the Majorana masses 
of sterile neutrinos [l^. Decays of a gauge singlet Higgs, 
playing also the role of the inflaton, were proposed as 
the origin of the relic population of sterile neutrinos in 



Ref. [6|]. The decays occur around the mass scale of 
the Higgs singlet, which, in that model, is taken to be 
in the sub-GeV range. The VEV of the singlet Higgs 
varies greatly from its mass. Here, however, we consider 
production of sterile neutrinos through decays of a sin- 
glet Higgs with mass and VEV both in the same energy 
scale [7, 8]. The requirement that a keV sterile neutrino 
constitutes all dark matter fixes this scale in the elec- 
troweak range, resulting in dark-matter particles with 
average momentum 3 times lower than in the Dodelson- 
Widrow production mechanism Q. In what follows, we 
will briefly review the definitions of free-streaming length 
and phase space density, and we will discuss their rel- 
evance to structure formation. We will then calculate 
these quantities for the suggested mechanism. The re- 
sults can be compared to observations. 



II. FREE-STREAMING AND PHASE SPACE 
DENSITY 

The free-streaming length Afs is the scale that corre- 
sponds to the cutoff of the power spectrum of the den- 
sity perturbations. In addition to suppressing the struc- 
ture on small scales, it may determine the filamentary 
structures in which the first stars form p^ . The free- 
streaming length is usually taken to be the comoving 
distance traveled by a particle within a Hubble time [31 ■ 
More detailed studies [19] have shown that Afs is the 
scale of the excitation mode that defines the crossover be- 
tween those perturbations that grow under gravitational 
instability and those that are Landau damped. Boy- 
anovsky ^20] calculated the marginal comoving wavevec- 
tor fcfs, for multicomponent dark matter, consisting of 
several species of collisionless particles that interact only 
gravitationally. At redshift z, k{s{z) is 
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where the summation extends over all species constitut- 
ing dark matter. In ([T|), all of the dependence on the 
redshift z is included in the prefactor (1 + z)^^. The 
mass density p and are meant to be evaluated at 

the present time. 

Here, we consider dark matter consisting of only one 
component. We define x = p/T to be the comov- 



2 



ing momentum. Then, in terms of the distribution 
function of the dark-matter population immedi- 
ately after the production has been completed, the free- 
streaming length, calculated at the present epoch (z = 
0), Afe = is 
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accounts for the dilution 



where the factor £ = — tt^? r 

and the redshift of the dark-matter particles, caused by 
the entropy release due to the decoupling of relativistic 
degrees of freedom since the time of production at T = 

-^prod ■ 

On the observational side, the Lyman-a forest power 
spectrum as measured by the Sloan Digital Sky Survey 
(SDSS) yields a limit of kfs ^ 2 Mpc"^ or Afs ^ 3 Mpc 
at redshifts z « 2 — 6, which translates into Afs ^1.1 Mpc 
today (z = 0) 

In order to study the phase-packing that leads to 
the appearance of central cores in galaxy halos, the 
phase space density, first introduced by Tremaine and 
Gunn [22| . is often employed. One has to discern be- 
tween the fine-grained phase density, which is the distri- 
bution function, and the coarse-grained one. Boyanovsky 
et al. generalized the coarse-grained phase density for 
arbitrary distribution functions of dark-matter particles. 
Here we will adopt, up to a change in units, this defini- 
tion: 



ally defined as Tc = ^ 27 a'^/AnG p. Then 
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where Voa — \f^o is the asymptotic velocity of the halo's 
rotational curve. 

Galaxy rotation curves provide a direct estimate of 
the enclosed mass density as a function of the distance 
from the galactic center, extending out to the flat por- 
tion. The latter yields an estimate for Wqo, or equiva- 
lently the dark- matter velocity dispersion a. Thus, ro- 
tation curves contain all the information needed to es- 
timate the phase space density of the core. Recent ob- 
servational data from the most dark matter dominated 
galaxies, the dwarf spheroidal galactic satellites, yield 

Qobs « 9 • 10-*^ - 2 • 10-4 ^ 

One may expect that at least some portion of dark 
matter retains the low entropy state, corresponding to 
the primordial phase density, and is naturally concen- 
trated in the center of the galaxy, forming its core. Sim- 
ulations show that this is the case for CDM and it can 
still be a valid approximation for WDM, provided that 
cores are formed without a lot of dynamical heating. The 
relation between the core radius and the velocity dis- 
persion cr, for a given dark-matter model, is a prediction 
of this assumption, which can be tested in halos. In- 
dependently of whether this approximation is good, the 
primordial phase density is bounded from below by the 
observable phase density, and this can be used to derive 
constraints for the mass of a dark-matter candidate. 

The phase space density ([3]) can be expressed in terms 
of the dark matter distribution function nix): 
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where N stands for the number density. V — N/ {p^)^^"^ is 
the dimensionless phase space density introduced in . 
For coUisionless particles and in the absence of self- 
gravity, Q (or V) is a Liouville invariant quantity p^ . 
Gravitational dynamics, however, can lead to decrease of 
the coarse-grained phase density. This signifies the en- 
tropy increase caused by dynamical heating during the 
non-linear gravitational clustering [3, IH, [2J, |2^ . 

The value of the primordial phase density Qq sets an 
upper limit to the dark-matter density p, for a given 
"velocity" dispersion, that is p ^ S^/^Qoct"^, where 

a = ■^fj2 ^ ^ ^ is the 1-dimensional velocity disper- 
sion in the nonrelativistic regime. For CDM, Qo — > oo, 
and there is no limit to the central density. For WDM, 
however, Qo is finite and halos form central cores rather 
than divergent cusps. Equivalently, this limit can be ex- 
pressed in terms of the core radius Tc, which, in the case 
of an isothermal sphere model for the halo, is convention- 
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III. STERILE NEUTRINOS PRODUCED FROM 
SINGLET HIGGS DECAYS 

Sterile neutrinos can be produced through oscillations 
of the active neutrinos Q . This scenario appears to be in 
conflict with a combination of the X-ray bounds [l^l and 
the Lyman-a bounds [2l| . It is possible to evade this con- 
straint, if the lepton asymmetry of the universe is greater 
than O(10~'^) or if decays of additional, heavier sterile 
neutrinos occur, introducing some additional entropy and 
contributing to cooling of dark matter p^ . More recent 
estimations of the free-streaming length and phase space 
density of sterile neutrinos produced via oscillations can 
be found in Refs. 

An alternative production mechanism arises if the 
Higgs sector is extended by a (real) gauge singlet field 
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S, which couples to the standard model Higgs, through 
a scalar potential, and to the right-handed neutrinos [S, 



- ViH,S) + h.c. 
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In this model, the Majorana masses of right-handed neu- 
trinos arise due to their Yukawa couplings to the S boson, 
through spontaneous symmetry breaking. Sterile neutri- 
nos can be produced from decays of the S bosons. This 
possibility has been explored in detail in Refs 0, [1| • The 
requirement that the resulting abundance of a keV ster- 
ile neutrino constitutes all of the observed dark matter 
fixes the energy scale of S bosons to be around the elec- 
troweak scale ~ 10^ GeV. Sterile neutrinos are produced 
with a nonthermal spectrum and with average momen- 
tum lower than the equilibrium value at the same temper- 
ature. Since their production occurs at temperatures just 
a factor of a few below the S boson mass, i.e. around the 
electroweak scale, sterile neutrinos are further redshifted 
by the entropy production due to the decoupling of the 
standard model degrees of freedom. We shall therefore 
call these sterile neutrinos "Chilled Dark Matter" . This 
chilling changes the quoted bounds for the sterile neu- 
trino mass, derived from small-scale structure formation 
considerations [2l|. In what follows, we compute the po- 
tentially interesting, for this purpose, quantities of free- 
streaming length and phase space density, discussed in 
the previous section. As in Ref. [sl, we distinguish be- 
tween two cases: in equilibrium and out-of-equilibrium 
singlet Higgs decay into sterile neutrinos. 



mass. Since the decay rate of the S bosons peaks at 
temperature Tprod = ms/fprod, with rpiod = 2.3 Q, the 
sterile neutrino population will be diluted by ^ ~ 33 [3], 
representing the decoupling of all of the standard model 
degrees of freedom. Then, the free-streaming length and 
phase space density, Eqs. and ([5]), are 
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B. Out-of-equilibrium decays 

If S bosons are weakly coupled to the standard model 
particles {a,uj ~ 0, Xhs ~ 10~^), they decouple early 
and decay at a later time, while out-of-equilibrium. In 
this case, the distribution function of the dark-matter 
sterile neutrinos produced is [§[: 
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where r = and ns{xs, r) is the S boson distribution 
function after freeze-out occurs, at rf = ms/Tf: 



ns{xs,r) = 



r -t- -\/a;| + 



X 



eV^%+^ „ 1 \rf + ^x% + J 



Here, A 



h-'Mp 



, where h'^ takes into account the decays 
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of S particles into all of the sterile neutrino species and 



A. Decays in equilibrium 

If the couplings in the scalar potential ([7]) of the Higgs 
singlet to the standard model Higgs are large enough, S 
bosons remain in thermal equilibrium down to low tem- 
peratures, and decay while in equilibrium. The distribu- 
tion function of dark-matter sterile neutrinos produced 
from decays of the singlet Higgs, in thermal equilibrium. 
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where / is the Yukawa coupling of the dark-matter 
sterile neutrino to the 5* boson (indices dropped 
for simplicity), ms is the S boson mass and 



10^^ GeV is the reduced Planck 



rprod 
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FIG. 1: The dimensionless time parameter r = corre- 
sponding to the peak of the sterile neutrino production rate 
from out-of-equilibrium S boson decays, vs A. For the range 
of A considered, S boson decays peak at temperatures only a 
factor of a few below its mass, that is before the decoupling 
of the QCD degrees of freedom. 



4 



\ 

\ \ ^ 






rr=2.5 

If = 1.5 



0.01 0.05 0.1 0.5 1 5 

A 




L . . 

0.01 0.05 0.1 0.5 1 5 

A 



FIG. 2: The free-streaming length and the phase space den- 
sity for sterile neutrinos produced from out-of-equilibrium S 

boson decays. is in units of (^)' ("^^) Mpc and 

is in units of (j^) (-^^)^ (km/s)-* • ^^18,11 A implies delayed 
decays, resulting in warmer dark matter. Sterile neutrinos 
produced through this mechanism constitute all of dark mat- 
ter for A ^ 0.1. 

the standard model fermions through the mass mixing 
with the standard model Higgs: 

a ^ 10^ GeV stands for the VEV of the 5* boson and 
B = is the branching ratio of S decays into the dark- 
matter sterile neutrino. 

The case of interest here is Xhs ~ 10~^. For lower 
values of Xhs, S bosons never come into equilibrium. 
This sets a lower limit on A, corresponding to the bb decay 
mode: A ^ 0.01 (where we used Xh,Xs < 1, imposed 
by the requirement of perturbativity of the potential). 
A determines how fast S particles decay. Small values 



of A lead to delayed decays and therefore warmer dark 
matter. Here, however, the minimum value of A ensures 
that S decays peak at temperatures just a factor of a few 
below its mass, as shown in Fig.[T] Since 5 bosons live in 
the electroweak scale, the dark matter will be produced 
before the decouphng of the QCD degrees of freedom, 
which implies ^ « 25 — 33. The free-streaming length 
and the phase space density of dark matter are shown in 
Fig. [2l For definiteness, we used ms = 200 GeV. Out-of- 
equilibrium S decays into a keV sterile neutrino produce 
the right amount of dark matter for A « 0.1 
IV. CONCLUSIONS 

We have calculated the free-streaming length and the 
phase space density of dark matter in the form of ster- 
ile neutrinos produced from decays of a gauge singlet 
in the extended Higgs sector. These two quantities may 
not capture the entire spectrum of density perturbations, 
but they serve as a useful set of parameters for comparing 
the model predictions with the data. We have considered 
two cases: in-equilibrium and out-of-equilibrium decays 
of the singlet Higgs. Comparison of the results, Eqs. ([9]) 
and pUj) with Fig. [21 shows that in-equilibrium decays 
result in colder dark matter than out-of-equilibrium de- 
cays. 

The quantities considered here have been calculated 
for other sterile neutrino production mechanisms in 
Refs. [l^llOl. Sterile neutrinos produced from the singlet 
Higgs decays have lower free-streaming length and higher 
primordial phase space density than neutrinos produced 
nonresonantly via the Dodelson-Widrow mechanism 
This is due to both their primordial nonthermal distri- 
bution and the chilling that occurs due to the entropy 
generation since their time of production. Compared to 
"cool" sterile neutrinos produced by net-lepton number 
driven resonant conversion [3], chilled sterile neutrinos 
have comparable free-streaming length and lower phase 
density. Chilled sterile neutrinos can account for all of 
dark matter, which would be in agreement with the cur- 
rent observational limits. The Higgs structure of this 
model can be probed at the Larg e Hadron Collider at 
CERN and at a linear collider [2g| . 

The author thanks D. Boyanovsky and A. Kusenko for 
helpful discussions. This work was supported in part by 
the DOE Grant No. DE-FG03-91ER40662 and by the 
NASA ATP Grant No. NAG 5-13399. 



[1] For review, see, e^g., A. Strumia and F. Vissani, 

afxrv:hep-ph /0606054 

[2] P. Minkowski, Phys. lett. B67 , 421 (1977); M. Gell- 
Mann, P. Ramond, and R. Slansky, Supergravity (P. van 
Nieuwenhuizen et al. eds.). North Holland, Amsterdam, 
1980; T. Yanagida, in Proceedings of the Workshop 



on the Unified Theory and the Baryon Number in the 
Universe (O. Sawada and A. Sugamoto, eds.), KEK, 
Tsukuba, Japan, 1979; S. L. Glashow, in Proceedings of 
the 1979 Cargese Summer Institute on Quarks and Lep- 
tons (M. Levy et al. eds.). Plenum Press, New York, 1980; 
R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 



5 



912 (1980). 

[3] S. Dodelson and L. M. Widrow, Phys. Rev. Lett. 72, 17 
(1994). 

[4] X. d. Shi and G. M. Fuller, Phys. Rev. Lett. 82, 2832 
(1999). 

[5] K. Abazajian, G. M. Fuller and M. Patel, Phys. Rev. 
D 64, 023501 (2001). A. D. Dolgov and S. H. Hansen, 
Astropart. Phys. 16, 339 (2002). T. Asaka, S. Blanchet 
and M. Shaposhnikov, Phys. Lett. B 631, 151 (2005); 
T. Asaka, M. Laine and M. Shaposhnikov, JHEP 0701, 
091 (2007); M. Shaposhnikov, Nucl. Phys. B 763, 49 
(2007); D. Gorbunov and M. Shaposhnikov, JHEP 0710, 
015 (2007); 

[6] M. Shaposhnikov and I. Tkachev, Phys. Lett. B 639, 414 
(2006). 

[7] A. Kusenko, Phys. Rev. Lett. 97, 241301 (2006). 
[8] K. Petraki and A. Kusenko, Phys. Rev. D 77, 065014 
(2008) 

[9] A. Kusenko and G. Segre, Phys. Lett. B 396, 197 
(1997); A. Kusenko and G. Segre, Phys. Rev. D 59, 
061302 (1999). M. Barkovich, J. C. D'Olivo and R. Mon- 
temayor, Phys. Rev. D 70, 043005 (2004); G. M. Fuller, 
A. Kusenko, I. Mocioiu, and S. Pascoli, Phys. Rev. D 
68, 103002 (2003); A. Kusenko, Int. J. Mod. Phys. D 13, 
2065 (2004). 

[10] E. K. Akhmedov, V. A. Rubakov and A. Y. Smirnov, 
Phys. Rev. Lett. 81, 1359 (1998); T. Asaka and M. Sha- 
poshnikov, Phys. Lett. B 620, 17 (2005). 

[11] P. L. Biermann and A. Kusenko, Phys. Rev. Lett. 96, 
091301 (2006); M. MapeUi, A. Ferrara and E. Pier- 
paoh, Mon. Not. Roy. Astron. Soc. 369, 1719 (2006); 
J. Stasielak, P. L. Biermann and A. Kusenko, Astro- 
phys. J. 654, 290 (2007); E. Ripamonti, M. Mapelh 
and A. Ferrara, Mon. Not. Roy. Astron. Soc. 375, 
1399 (2007); J. Stasielak, P. L. Biermann, A. Kusenko, 
arXiv;astro-ph/0701585; Acta Phys. Polon. B 38, 3869 
(2007) 

[12] T. Asaka, A. Kusenko and M. Shaposhnikov, Phys. Lett. 
B 638, 401 (2006). 

[13] G. Kauffmann, S. D. M. White and B. Guiderdoni, Mon. 
Not. Roy. Astron. Soc. 264, 201 (1993); A. A. Klypin, 
A. V. Kravtsov, O. Valenzuela and F. Prada, Astro- 
phys. J. 522, 82 (1999); B. Moore, S. Ghigna, F. Gov- 
ernato, G. Lake, T. Quinn, J. Stadel and P. Tozzi, ApJ 
524, L19 (1999). B. Willman, F. Governato, J. Wad- 
sley and T. Quinn, MNRAS, 355, 159 (2004); P. Bode, 
J. P. Ostriker and N. Turok, Astrophys. J. 556, 93 (2001). 
P. J. E. Peebles, ApJ, 557, 495 (2001); A. R. Zent- 
ner and J. S. Bullock, Phys. Rev. D 66, 043003 (2002); 
J. D. Simon, A. D. Bolatto, A. Leroy and L. Blitz, As- 
trophys. J. 596, 957 (2003); F. Governato et al, As- 
trophys. J. 607, 688 (2004); G. Gentile, P. Salucci, 
U. Klein, D. Vergani and P. Kalberla, Mon. Not. Roy. 
Astron. Soc. 351, 903 (2004); J. Kormendy, M. E. Cor- 
nell, D. L. Block, J. H. Knapen and E. L. AUard, As- 
trophys. J. 642, 765 (2006); M. I. Wilkinson et al, 
arXiv;astro-ph/0602186; L. E. Strigari, J. S. Bullock, 
M. Kaplinghat, A. V. Kravtsov, O. Y. Gnedin, K. Abaza- 
jian and A. A. Klypin, Astrophys. J. 652, 306 (2006); 

[14] J. J. Dalcanton and C. J. Hogan, Astrophys. J. 561, 35 
(2001); 

[15] D. Boyanovsky, H. J. de Vega and N. Sanchez, Phys. Rev. 

D 77, 043518 (2008) 
[16] Y. Chikashige, G. Gelmini, R. D. Peccei and M. Ron- 



cadeUi, Phys. Lett. B 94, 499 (1980). Y. Chikashige, 
R. N. Mohapatra and R. D. Peccei, Phys. Lett. B 98, 
265 (1981). 

[17] L. Gao and T. Theuns, Science 317, 1527 (2007). 
[18] E. W. Kolb, M. S. Turner, The Early Universe , Addison- 
Wesley, (1990). 

[19] I. H. Gilbert, Astrophys. J. 144, 233 (1966); ibid, 152, 
1043 (1968); J. R. Bond, A. S. Szalay, Astrophys. J. 274, 
443 (1983). 

[20] D. Boyanovsky, Phys. Rev. D 77, 023528 (2008) 
[21] M. Viel, et al., Phys. Rev. D 71, 063534 (2005); U. Seljak, 
A. Makarov, P. McDonald and H. Trac, Phys. Rev. Lett. 
97, 191303 (2006); M. Viel, et al., Phys. Rev. Lett. 97, 
071301 (2006); A. Palazzo, D. Cumberbatch, A. Slosar 
and J. Silk, Phys. Rev. D 76, 103511 (2007); M. Viel, 
G. D. Becker, J. S. Bolton, M. G. Haehnelt, M. Ranch 
and W. L. W. Sargent, Phys. Rev. Lett. 100, 041304 
(2008) 

[22] S. Tremaine and J. E. Gunn, Phys. Rev. Lett. 42, 407 
(1979). 

[23] D. Lynden-Bell, Mon. Not. Roy. Astron. Soc. 136, 101 
(1967); S. Tremaine, M. Henon, D. Lynden-Bell, Mon. 
Not. Roy. Astron. Soc. 219, 285 (1986); 

[24] C. J. Hogan and J. J. Dalcanton, Phys. Rev. D 62, 063511 
(2000). 

[25] S. Peirani et. al, Mon. Not.R. Astron. Soc. 367, 
1011 ( 2006); S. Peirani, J. A. de Freitas Pacheco, 
'astro-ph/0701292 ; E. Romano-Diaz et.al, Astrophys. J. 
637, L93 (2006); E. Romano-Diaz et.al, Astrophys. J. 
657, 56 (2007); Y. HofTman et.al arXiv:0706.0006 

[26] L. E. Strigari, J. S. Bullock, M. Kaphnghat, 
A. V. Kravtsov, O. Y. Gnedin, K. Abazajian and 
A. A. Klypin, Astrophys. J. 652, 306 (2006); L. E. Stri- 
gari, M. Kaplinghat and J. S. Bullock, Phys. Rev. 
D 75, 061303 (2007); G. Gilmore, M. I. Wilkinson, 
R. F. G. Wyse, J. T. Kleyna, A. Koch, N. W. Evans 
and E. K. Grebel, Astrophys. J. 663, 948 (2007); 
R. F. G. Wyse and G. Gilmore, arXiv:0708.1492 

[27] K. Abazajian, G. M. Fuller and W. H. Tucker, Astro- 
phys. J. 562, 593 (2001); A. Boyarsky, A. Neronov, 
O. Ruchayskiy and M. Shaposhnikov, Mon. Not. Roy. 
Astron. Soc. 370, 213 (2006); A. Boyarsky, A. Neronov, 
O. Ruchayskiy and M. Shaposhnikov, JETP Lett. 83, 
133 (2006); A. Boyarsky, A. Neronov, O. Ruchayskiy, 
M. Shaposhnikov and I. Tkachev, Phys. Rev. Lett. 97, 
261302 (2006); S. Riemer-Sorensen, S. H. Hansen and 
K. Pedersen, Astrophys. J. 644, L33 (2006); K. Abaza- 
jian and S. M. Koushiappas, Phys. Rev. D 74, 023527 
(2006); C. R. Watson, J. F. Beacom, H. Yuksel 
and T. P. Walker, Phys. Rev. D 74, 033009 (2006); 
K. N. Abazajian, M. Markevitch, S. M. Koushiappas and 
R. C. Hickox, Phys. Rev. D 75, 063511 (2007); A. Bo- 
yarsky, J. Nevalainen and O. Ruchayskiy, Astron. As- 
trophys. 471, 51 (2007); A. Boyarsky, O. Ruchayskiy 
and M. Markevitch, arXiv;astro-ph/0611168 , S. Riemer- 
Sorensen, K. Pedersen, S. H. Hansen and H. Dahle, Phys. 
Rev. D 76, 043524 (2007); A. Boyarsky, J. W. den 
Herder, A. Neronov and O. Ruchayskiy, Astropart. 
Phys. 28, 303 (2007 ); H. Yuksel, J. F. Beacom and 

C. R. Watson, lajXi7 :0706.4084 [astro-ph]; A. Boyarsky, 

D. lakubovskyi, O. Ruchayskiy and V. Savchenko, 
arXiv:0709.2301 [astro-ph]. 

[28] T. Binoth and J. J. van der Bij, Z. Phys. C 75, 17 (1997); 
A. Datta et al., Z. Phys. C 72, 449 (1996). H. Davoudiasl, 



6 



T. Han and H. E. Logan, Phys. Rev. D 71, 115007 (2005); 
M. J. Strassler and K. M. Zurek, arXiv:hep-ph/0605193 
D. O'Connell, M. J. Ramsey-Musolf and M. B. Wise, 
Phys. Rev. D 75, 037701 (2007); V. Barger, P. Langacker 
and G. Shaughnessy, Phys. Rev. D 75, 055013 (2007); 



S. Profumo, M. J. Ramsey-Musolf and G. Shaughnessy, 
JHEP 0708, 010 (2007); V. Barger, P. Langacker, M. Mc- 
Caskey, M. J. Ramsey-Musolf and G. Shaughnessy, Phys. 
Rev. D 77, 035005 (2008) 



